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Manufacturing Cost and Strength Reliability Effects 
on the Optimal Design of Hand Prosthesis 

Mechanisms 

Emmanuel Ngale Haulin, Nestor Tsamo 


Abstract — The middle finger mechanism of a hand prosthesis 
is optimally designed with respect to five design criteria: motion 
posture error, mechanical error, maximum driving torque, 
strength reliability and manufacturing cost of the mechanism. It 
is shown that the optimal design variables obtained are different 
from one design criterion to another when optimized 
individually. Otherwise, the introduction of the global 
manufacturing cost to the multiobjective optimization of the 
first three criteria on the one hand and of the first four criteria 
on the other hand, gave the same optimum variables obtained 
when minimising simultaneously only the three first criteria. 
Therefore, the global manufacturing cost has a real impact on 
the whole process of the mechanism optimization due to the 
stochastic nature of the trajectory of the cutting tools during the 
whole process of manufacturing. 

Index Terms — hand prosthesis, manufacturing cost, 
mechanisms, optimal design, optimization criteria, strength 
reliability. 


I. Introduction 

Assuming that all fingers of a natural hand, except the 
thumb, are functionally and mechanically similar during the 
generation of the various modes of a basic gripping of the 
hand, a simple mechanism with one degree of freedom is often 
used [1,2,3,4] in order to proceed to a multiobjective 
optimization of the middle finger mechanism only. 

Among many mechanisms and criteria used by many 
authors to optimize the middle finger, we choose the across 
four-bar mechanism usually optimized by the means of 
motion posture error, maximum driving torque [1,2,3,4] and 
mechanical error due to dimensional tolerances [1,2] which 
have a significant effect to that mechanism. In order to obtain 
a realistic optimal design, two additional criteria, namely 
strength reliability and manufacturing cost, will be also used 
and their effect in the optimization process will be studied. 

In a first approach, each criterion will be optimized 
individually. Then, the three first criteria will be optimized 
simultaneously. Strength reliability and manufacturing cost 
will be added individually and finally, we will carry out the 
optimization of all the five criteria simultaneously. 

II. MIDLE FINGER DESCRIPTION 

A. Finger model 

The type of movement achieved by a natural hand is a 
balance-balance one. Therefore, the planar one degree of 
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freedom and across four-bar mechanism used to model the 
finger of the prosthesis is also a balance-balance mechanism 
(fig. 1 [2,5]). It is equipped with the following elements: 

• a driving system consisting of a cable winded around a 

pulley. 

• a fixed bar QM on the palm of hand; 

• a driving bar AM activated by the driving system and 

allowing phalanges to bend by the means of the driven 
bar QB and the junction bar AB; 

• a return spring allowing not only to maintain a minimum of 

tension in the cable, but also to make sure of the finger’s 
return to its opening or operating position when the cable 
is unwinded. 


Interphalangeal 



Figure 1: Across four-bar mechanism of the middle finger 
B. Description of the middle finger's mechanism 


Fig. 2 shows the middle finger’s parameters related to 
the across four-bar mechanism equipped with its driving 
system. 



P3 ’: effective length of the third phalanx P3; 

EXC: eccentricity of the articulated joint between the 
first and second phalanges; 
re length of driven bar QB related to first phalanx; 
r 2 : length of junction bar AB related to second phalanx; 
r 3 , r 4 : respectively lengths of driving bar AM and fixed 
bar QM; 
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TP1, cp, \|/: bending angles of phalanges PI, P2 and P3 
with respect to horizontal palm plan; 

TP2, TP3: bending angles of phalanges P2 and P3 with 
respect to PI and P2; 

0i, 0 2 , 03 , 04 : respectively angles of bars QB, AB, AM 
and QM with respect to palm plan; 

0 A : angle between junction bar AB and second phalanx 
axis; 

x T , y T : connection T coordinates with respect to M; 
x c , Yc- pulley center coordinates with respect to Q; 
d: pulley diameter; 

p: traction angle; y: position of contact point S between 
pulley and cable; 

P: applied force to the finger end with respect to the third 
phalanx P3; 

F t : tension in the cable; 

M m : necessary driving torque applied to the driving bar 
AM; 


C. Mechanism synthesis 

At any mechanism position i, Freudenstein’s relation 
[1,2,4,5] between 0n and 0 2 i is : 

k\ cos(<9 i; -0 4 )+k 2 cos(<9 2 , -0 4 )-k 3 = cos (0 U -0 2i ) (1) 

2,2 2,2 

— —; k 2 = ~ — — + r *~ r 3 +r * 

r 0 r, 


with 


2r \h 


r x =4p\ 2 +EXC 2 
0 U =TPl i -arctg(- 


\e 2l =0 A+ TP\ l+ TP2, 

Assuming that r 1? 0 A , 0 4 , 7 positions TPli and TP2i are known, 
the mechanism synthesis is then performed using the least 
square method [1,2,4]. 


III. Optimization of the middle finger’s mechanism 


A. Mechanism parameters 

The force P applied at the end of middle finger’s is equal 
to 45N and makes an angle 0 P =9O° with the phalanx P3 axis. 
The maximum tension developed in the cable is F Tmax =400N. 
Table 1 gives overall dimensions of middle finger’s 
mechanism with DPI, DP2 and WR3 which are respectively 
diameters of phalanges PI, P2 and width of links. 

Table 1: Linear and angular dimensions of middle finger’s 
mechanism [1,2] 


(mm) 

DPI 

DP2 

WR3 

PI 

P2 

P3’ 

EXC 



15 

13 

5 

43.5 

29 

13.5 

2 


o° 

TPf 

TP1 2 

TP1 3 

TP1 4 

TP1 5 

tpi 6 

TP1 7 



18 

29 

40 

51 

61 

73 

84 


o° 

TP2i 

TP1 2 

TP2 3 

TP2 4 

TP2 5 

TP2 6 

TP2 7 

TP3 


14.274 

18.027 

24.022 

32.142 

41.416 

55.649 

75.670 

30 


B. Design variables 

From the mechanism parameters and with the aim to 
introduce the driving system in the optimization process, 
vector X of the independent design variables is: 
X=[x 1 ,x 2 ,x 3 ,x 4 ,x 5 ,x 6 ,x 7 ,x 8 ,x 9 ] 

= [0 A ,r 2 , r 3 ,r 4 ,0 4 ,t m ,x T ,y T ,d\ (2) 


with t m = thickness of different bars and 
e A = 0 2 . -{TP2 l + TP ],) . Table 2 gives design variables 
bounds. 

Table 2: Design variables bounds related to middle finger 
mechanism [1,2] 


(mm) 

r 2 

r 4 

e 4 (°) 

tm 

x T 

y T 

d 

min 

5 

5 

360 

0.794 

10 

-5 

10 

max 

15 

7.5 

380 

1.588 

20 

5 

17 


C. Objective functions 

Five performance criteria are considered in the process 
of the mechanism optimization: 

• motion posture error or 

structural error on the bending angle of the second 
phalange that would be minimized to ensure that the 
prosthetic finger should pass by the prescribed positions 
[1,2,4]; 

• mechanical error on the 

bending angle of the second phalange, due to 
dimensional tolerances and clearances on the 

articulations, that would be minimized[1,2]; 

• maximum shaft driving 

torque that would be minimized and able to 

counter-balance a grasping force applied at the end of the 
middle finger [1,2,4]; 

• strength reliability of the 

mechanism that would be maximized to ensure the 

mechanism strength during grasping, holding or pinching 
operations; 

• manufacturing cost of the 

mechanism that would be minimized to evaluate its effect 
on the design variables; 

Objective function related to the structural error on bending 
angle TP2 

.f\(x)= X -Yid 2lc -d 2ld ) (3) 

1=1 

where 0 2 id = 0 A + TPf + TP2i = Xj + TPf + TP2i is the 
desired angle and 0 2ic the computed or real angle. 

Objective function related to the mechanical error on 
bending angle TP2 

/2 (,)=i£(A<y 2 (4) 

1=1 


with A 0 2i = 3<j0 where standard deviation is : 



and cj = 0.001a, 

u j J 


Objective function related to the maximum shaft driving 
torque 
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f 3 (x)=max(M P ) 


(5) 


with 


M P = 


M m d 


x 9 M m 


2{x t sin/? + y T cos/?) 2{x 7 sin/? + x 8 cos/?) 


M P is the driving torque applied to the shaft while M m is 
torque reduced to point M of the driving bar [1,3,6], assuming 
that frictional forces at joints are negligible. 

Objective function related to the strength reliability of the 
mechanism 

Since the failure of any link causes the breakdown of the 
whole mechanism, the across-four bar mechanism is 
comparable to a system made up of four components 
assembled in a serial configuration. Therefore, its strength 
reliability is the product of the reliability of driven, junction 
and driving bars R b R 2 , R 3 [ 6 , 8 ] given by the mathematical 
expression: 

f 4 (x) = - Ri.R 3 .R 3 ( 6 ) The minus sign is used to 

minimize the objective function f 4 (x) in order to maximize the 
strength reliability of the mechanism. The strength reliability 


of any element is R = 



/ \ 

7.8125 

1- 

1.29 


\u\ 




v 1 1 / 



[ 8 ]. Assuming that the 


maximum stress p induced in a bar and the strength S of 
material used are lognormally distributed, the statistical 


parameter is u - 


In /?-ln S 


* 


[ 1 , 2 , 8 ] where 


In p 


+ CT, 


q = tgd x ; c 2 =p x = -P cos(o P + if/}, 

C 3 = P y = Psm(0 P +y/\ C 4 = r 2 sm0 2 ; 

with ]c 5 =r 2 COS0 2 ; C 7 = C, Q \ 

C 6 = P y (P3'cosy/ + P2cos(p-EXCsm(p)+ 
P X (P3' sin y/ + P2sin (p + EXCcostp) 




Figure 3: Free body diagram of junction bar AB 


Driven bar QB (fig. 4) 

t 


A - r maxl - 


In S 


Igv 

2 


Fbx + Fby 
2t yy ,DPl 


( 8 ) 


\na = \na-0.5(j] 


In a"> u ln a 


= In 




+1 


cr =0.001 


- 1 

da ^ 

N 

_l 

1 ^ 


+1 

7=1 


da 

dx T 




a is the mean value of a, a a is the standard deviation of stress 
or strength a. Using the maximum shear theory, p = x max is the 
maximum shear stress induced in element and S = S y /2 where 
S y = 1345MPa is the yield strength of material used. The 
maximum shear stress p induced in each element is 
determined from the free body diagram of each bar: 

Junction bar AB (fig. 3) 


Pi 


with 


L max 2 


2 


1 


2t„DP2 


6 M 


f max | 


DP 2 


(7) 


N = -P[cosO P cos(y/ -(p)-sm 0 P sm(y/-(p)\ 


M 


f max 


= -P\P 3 sin 6 P + 


C? 


Fbx ~ ~ ’ Fby ~ Ci ’ Fax - Fi 

c i 


P 2 (cosO P sin(^ - cp)+ sin 0 P cos(i// - <^))] 
bx +C 2 \ F ay = E by - C 3 


i F Qy 


rQx 





Figure 4: Free body diagram of driven bar QB 

Driving bar AM (fig. 5) 

W Fmx' - Fax' 


P 3 - r max3 - 


2 


1 


2t m WR3 

1 


6 \M 


/max 


AX' 


WR3 
6 | M 


f max | 


2x 6 WR3 


if f mx'— r 


AX' 


WR3 


AX' 


-1 


2t m WR3 

-1 


Fmx' + 


6 M 


/max | 


2x 6 WR3 


WR3 
6 | M 


f max 


[ MX' 


WR3 


if (y T < 0 and F MY < < 0) or (y T > 0 and E MT > 0) 


\M 


—\Xt\Emv' — x 7 \F 


f max _ | A 7T MY 


MY' 


(9) 
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if (y T < 0 and F ur > 0) or ( y T > 0 and F MY . < 0) 
i i wi 

|^/max| = + C0S /? = X 1 \^MY' | + | X 8 \^T C0S /? 


th 


p _ r 3 Fay' ___ r 3 Fay' _. 

T (x T sm/3 + y T cosfl) (x 7 sin/? + x 8 cos/?)’ 

Fmy' ~ -F t sin fi-F Ar 

Fmx' = Fj cos/]-F ax F AX . = F AX cos0 3 + F AY sin 0 3 ; 
Fay' ~ Fay cos6 3 - F ax sin 0 3 



Figure 5: Free body diagram of driving bar AM 


Objective function related to the global manufacturing cost 
The assembly drawing of the middle finger (fig. 6) is done 
according to parameters of the mechanism described in 
section II.B and taking into account the width WR3 of the 
driving link, the diameters DPI of the first phalanx and DP2 
of the second phalanx. 



Figure 6: Assembly drawing of the middle finger with its 
across four-link mechanism 


1. Palmer frame 

2. Guiding screw 01 

3. Guiding screw 02 

4. Involved bar 

5. Guiding screw03 


6. Junction bar 

7. Driving bar 

8. Pressure screw 
9. Cable 

10. Pulley 


Basic assumptions 

• Stainless steel pieces have overall fixed dimensions 2x5xL 

and 2x11.5xL where L is the length; 

• Unit manufacturing cost is expressed in terms of 

manufacturing time and machine tool hourly cost; 

• Because of the programmable machine tool used, machine 

hourly cost is constant; 

• Unit manufacturing costs of machining operations such as 

turning, drilling, boring, facing and undercutting are 
equal; 


• Each manufacturing dimension is obtained in only one final 

operation; 

• Machining linear velocity V is in m/min; 

• Tools store room is equipped with necessary cutting tools 

used on the programmable machine tool; 

• Thermal deformations of machine-piece-tool set are 

stabilized. 

Global manufacturing cost 

The present economic constraints impose to the enterprises of 
production, the maximal profit that is the production at a 
minimal cost and at a minimal time [11,12,13]. The advanced 
methods of machining as the use of the manufacturing 
Computer Aided Design, improved the productivity of the 
machining operation meaningfully. The prostheses are 
manufactured in very small set, therefore, the manufacturing 
Computer Aided Design is appropriate for this purpose [14]. 
Let’s consider : 

C g : The global manufacturing cost of across four-bar 
mechanism ; 

C(i): The global manufacturing cost of the i th bar of the 
structure of the mechanism, where i=l,2,3 which 
corresponds respectively to the involved bar, driving bar 
and junction bar ; 

C u (ij): The manufacturing unitary cost of the j th Manu6 
facturing operation on the i th bar of the mechanism, where 
j=l,2,3,4,5 corresponding respectively to facing milling, 
lateral milling, drilling, counter boring and tapping ; 

D ( ij) : Cutting tools used relatively at i th bar and j th 
manufacturing operation (Milling cutter 2 cuttings 020, 
cutter 2 cuttings 06, Drill 02, Drill 04, piloting Milling 
cutter to counter boring 02x6 and Machining tapping 
M2); 

L(ij) — ^s(ij) - L C (ij) Lp(ij) L^/(ij) L t(i j) . Tools trajectories 
relatively to the facing milling, lateral milling, drilling, 
counter boring and tapping of the the i th bar of the 

mechanism and the j th manufacturing operation; 
OP(i)=2(Xo(i) + Yo(i))=XY(i): Starting machine programme = 

Starting piece programme of the i th bar of the mechanism ; 

V, N, f, k, C : respectively linear speed, mass-production 
set, advance in turn/tooth, Taylor constant function of 
material of the cutting tool, Taylor constant function of 
material of the working piece (NF E 66-505) ; 

Ti T 2 T 3 : Holes 1,2, and 3 ; 

A, F, R, L, E, S : respectively technical amortization, 
financial expenses, maintenance expenses and of repair, 
expenses of local or of clutter, expenses of 
energy and wage costs and social. 

C m : hourly machine cost (Advanced Manufacturing 
Computer Aided Design) 

C s : Cost cutting tools ; 

C p? C T , C Tm , C 0 : respectively sum of the expenses related 
directly to the preparation, cutting duration, out of cutting 
time and the cutting tool; 

Ap : global rate of exploitation of the preparation section, 
general expenses and labour ; 

Expression of the unit cost of machining (C u(i j)) 

The production to the minimal cost and to the maximal profit 
called on the mathematical models translating the laws of 
wear of the tools according to the cutting times also named 
simplified Taylor’s law [15]. The manufacturing unitary cost 
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of the j th manufacturing operation on the i th bar of the 
mechanism is C u ^ — C p + C T + C Tm + C 0 with 

c = L X ^. C gm .. %) x,rx % . 

p N 60 ’ T 60 10 3 xVx/ 


C - y 
Tm 60 




c„=c. 


L(ij) x k x A, 


10 3 x f xC~ K 


My-tf-i. 


TV 


C m — A+F+R+L+E+S 

The final expression of the cost of a machining pass (turning, 
drilling, milling...) is: 

C u(a) = C F +7iL [ij) D [ij) M 


with 


M = 


CL 


6 xl0 4 xVx / 


\ f 

+ 




10 3 x / x C~ K 


V 


-K-l 


C F - Cp+ C Tm N 


Manufacturing cost of the junction bar 
Fig. 7 shows the manufacturing drawing of junction bar. It is 
manufactured in four basic operations : facing milling L s(2 i), 
lateral milling L c(2 2 ), drilling L p(23) , and counter boring L/ (24) . 
The unit manufacturing cost related to each of these 
operations are : 

4( 3 l) = 4 + ^(31)^(3 1 )^’ 4( 32 ) = 4 + L c ( 32)^(32)^’ 
4(33) = Cp + 4( 33)D(33)M’ 4( 34) = C F +L / ( 34 )Z)( 34 )M; 

Therefore, the manufacturing cost of the junction bar is : 

C 3 = 4 C F + 44 (31 )/) (31 ) + 4( 32 )£>( 32 ) + 4(33)433) + 4(34)434) V 


Cutting tool profile 



h 

■SL 





Figure 7: Junction bar 


Manufacturing cost of the driving bar 

Fig. 8 shows the manufacturing drawing of driving bar. The 
trajectories of the cutting tools are established according to 
the working drawings [7]. This bar 1 is manufactured in four 
basic operations : facing milling L s( n), lateral milling L c(12) , 
grilling L p(13) , counter boring L/ (14) . The unit manufacturing 
cost related to each of these operations are : 

4(n) = 4 + ^(ii)4n)^; 4(i2) = 4 + 4 (i 2)4 i2 )^ ; 

4( 13 ) = 4 +4(i3)4i3) M ; 4( 14 ) = 4 +4 i4 )4 i4 )^ ; 

erefore, the manufacturing cost of the driving bar is : 

Q = 4 C F + 44(11)0(11) + 4( 12 )D(12) + 4(13)433) + 4(14)414)]^ 



Manufacturing cost of the driven bar 

Fig. 9 shows the manufacturing drawing of driven bar. It is 
manufactured in five basic operations : facing milling L s(21) , 
lateral milling L c(22) , drilling L p(23) , counter boring L/ (24) , 
tapping L t(25) . The unit manufacturing cost related to each of 
these operations are: 


u(2l) =4 + ^4(2l)4 2 l) M ’ 4(22) =4 + L c(22 ) D (22) M \ 
u( 23 ) = 4 +4(23)4 23 ) M ; 4(24) = 4 + 4 24 )4 24 )^’ 


4(25) ~~ 4 +425)4 2 5)^ ; 

Therefore, the manufacturing cost of the driven bar is: 


C9 — 5Cr + 7T 


4(2l)4 2 l) + 4 22 ) 42 ) 

+ 4 ( 23 ) 433 ) + 4 ( 24 ) 424 ) + 4(25 ) 4 5 ), 


M 



Global manufacturing cost of the mechanism Cg 
The global manufacturing cost of the mechanism is 
= Cj + C 2 + C 3 = 13C^ + tM x LD (10) 

With 


LD = 


4( 1 1)4 1 ) + 4 12 )4 12 ) + 4(13)4 13 ) + 4(l4)4 14 ) a 
+ ^(2l)4 21 ) + L c{ 22 ) D ( 22 ) + L p(23) D (23) + 4 24 4 24 ) 

+ L l(25) D (25) + 4 ( 3 l ) 4 3 l ) + L c(32) D {32) + L p{33) D {33) 

+ L l{34) D (34) j 

is in US DOLLARS and the dimensions of all mechanism 
pieces are in millimeter. 


While taking into account the dimensional tolerances Ar ; . to 
the bars, clearances at the articulations and an 

addimensionnal variable k [2] : 

Ar, ARc, ~ 

-^- =-— = = O.OOlr. => Ar.- = 3.10 r.- (11) T 

3 3 J J J J 

ables 3, 4a and 4b give respectively constant expressions and 
bars profile coordinates used in the global manufacturing cost 
Cg. 
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Table 3 : Constant expressions 


Tools Trajectories 

Constants 

Lsii = Qi+2x 3 

Qi=4(ei+Ri) 

L c i 2 = Q 2 + x 3 

Q 2 =27iRi + 1 2 + 1 3 +Ci+ c 2 + c 3 

Lpi 3 - Q 4 + Zi+Z 3 

Q 3 =OPi + 3 81 

L/i 4 = Q3+2(x 7 +x 8 )+3x 6 

Q 4 =OP 2 + Zi+Z 3 

Q 5 -OP 2 + 4 s 2 + 2 R 2 + 1 4 + I 5 

Ls21=Q5+ r l 

+ Z s + 1 6 

Lc 22 = Q 6 + r l 

Qg— OP 2 +2s 3 +tt R 3 + 7i R 4 

Lp 23 = Qll + Ql2+2Qi 3 

+ 1 9 + C 5 + 1 8 + C 4 + lio 

+2Zp +x 4 

Q 7 — OP 2 + 2s 3 +2 s 4 +2R 2 

+2 (rj+x 6 ) with 

Qs= OP 2 + 2 s 3 

x 4 =X 4 [cos (271- G 4 ) 

+2R 2 +2(s 4 +Z 2 ) 

-\-sin(2n-G 4 )] and 

Q 9 — OP 2 + 2s 2 +2R 2 +28 4 

0 4 — x 5 

Q 10 — OP 2 + 28 2 +2 s 4 +2Z 4 

L/2 4 =27+2s+229+2Z/? 

Qn= OP 2 + 28 3 +2 s 4 +2R 2 

+ x 4 +2 (rj+x 6 ) 

Q 12 — OP 2 +2s 3 + 2R 2 

L t 25 = Qio+2Z 4 

+2(s 4 + Z 2 ) 

Qi 3 — OP 2 + 28 2 +2R 2 +2s 4 

L s3 i= Qi 4 

Qi 4 = OP 3 +2(2s 3 +R 4 +P 2 + P 3 ) 

L c3 2= Ql5 

Q15= OP 3 +2s 4 +7iR 4 +P 2 

Lp 33 =Ql 6 +Ql 7 + 2 (x 6 +X 2 ) 

+ C 7 + C 8 + C 9 + Cio"l" Cn 

+ 2 x 6 

+ Ci 2 + li 3 

L /34 =2x 6 +Q 18 

Qi6— OP 3 + 28 3 + 2R 4 + 2s 4 

Qi 7 - OP 3 + 28 3 + 2R 4 + 2s 4 

Qi 8 = OP 3 + 28 3 + 2R 4 + 2s 4 


Table 4a : Bars profile coordinates 


Bars 

Opera¬ 

tions 

relatives profile coordinates 

Driving 

Facing 

milling 

Lateral 

milling 

Drilling 

Counter 

boring 

X 01 

Y 01 

ii 

I 2 

1 3 

Cl 

C2 

C 3 


30,2 

35,2 

x 3 

29,6 

29,6 

6,7 

22,4 

6,7 


£ i 

tm 

Xt 

y T 

r 3 

X 10 


Ri 


3 

x 6 

X 7 

x 8 

x 3 

19,2 


11 


Involved 

Facing 

milling 

Lateral 

milling 

Drilling 

Counter 

boring 

Tipping 

Zi 

z 2 

r 3 

£ 2 

r 2 

u 

i 5 

k 

R s 

1 

1,2 

X 3 

5 

11 

62,1 

5,6 

5,6 

7 

r 3 

r 4 

£ 3 

c 4 

c 5 

18 

I 9 

ho 


10 

4,2 

5,6 

14 

0,8 

62,1 

62,1 

0,8 


x 02 

Y 02 

Z 2 

z p4 

Z 4 

x 4 

x 6 



23,2 

22,8 

2 

7 

1 

x 4 

x 6 




Table 4b : Bars profile coordinates 


Bars 

Opera¬ 

tions 

relatives profile coordinates 

Junction 

Facing 

milling 

Lateral 

milling 

Drilling 

Counter 

boring 

X 03 

Y 03 

84 

r 4 

cl 

c7 

c 8 

c9 

clO 

cl 1 

5,7 

5,7 

3 

9,9 

7,5 

8,5 

3,8 

16,8 

35,5 

10,9 












Objective functions related to multiobjective optimization 

f 6 (x)= [fi(x),f 2 (x),f 3 (x)]; 

f?(x) = [fi(x), f 2 (x), f 3 (x),f 4 (x)]; 

f 8 (x) = [fi(x), f 2 (x), f 3 (x),f 5 (x)J ; 

f 9 (x) = [fi(x), f 2 (x), f 3 (x),f 4 (x), f 5 (x)]. 

Each of these multiobjective functions consists of a set of 
objective functions which are simultaneously optimized. 

D. Constraint functions 

These functions, expressed in terms of design variables 
are equality constraints hi(x) = 0 and inequality constraints 
gi(x)< 0 according to Matlab optimization toolbox R2015a. 

Equality constraints 

Using Freudenstein’s relation and least square method, three 
equality constraints [1,2,4,5] related to mechanism synthesis 
with reference to seven positions given in table 1 are obtained. 
Otherwise, one additional equality constraint is introduced by 
the means of security factor related to the design of the driving 
bar. 


Inequality constraints 

Eighteen inequality constraints [1,2] used are respectively 
related to transmission angle, mechanism dead points, 
restriction of driving link to be within the phalanges PI and 
P2, maximum value of cable tension, maximum driving shaft 
rotation, connection cable T keeping inside phalanx PI during 
a cycle of mechanism’s closure, cable keeping within palm 
and finger. 

In addition to design bounds given in table 2, angle 0 2 i must to 
be included between 20° and 115°[1,2]; therefore 0 A is 
bounded by -12.274° and 82.726°. 

IV. RESULTS AND DISCUSSION 

A. Optimization results and interpretation 

From three different starting points, each of objective 
functions fi(x), f 2 (x), f 3 (x), f 4 (x) and f 5 (x), subjected to four 
equality and eighteen inequality constraints, was optimized, 
taking into account design variables bounds, by the means of a 
nonlinear programming algorithm developed in Matlab 
optimization toolbox R2015a. 

Using the optimum values of each single objective function, 
the multiobjective optimization of f6(x), f 7 (x), f 8 (x), and f 9 (x) 
subjected to the same constraints was carried out by the means 
of Matlab goal attainment method [1,2] developed in Matlab 
optimization toolbox R2015a. 
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Table 5: Optimum design variables and functions 



fi(°> 2 

f 2 (°) 2 

f 3 

(Nmm) 

u 

f 5 

(US $) 

f« 


fs 

u 


0.1595 

13.6892 

1307.978 

0.9999 

494.845 





Xi(°) 

82.726 

73.409 

72.257 

82.518 

82.726 

78.193 

82.152 

82.063 

82.063 

x 2 (mm) 

5.000 

5.148 

5.000 

5.065 

5.604 

5.060 

5.094 

5.060 

5.060 

x 3 (mm) 

37.153 

36.741 

36.954 

37.046 

36.232 

37.035 

36.987 

37.035 

37.035 

x 4 (mm) 

5.856 

7.003 

6.872 

5.973 

6.724 

6.019 

6.057 

6.019 

6.019 

X 5 (°) 

375.063 

380.000 

380.000 

375.415 

377.091 

375.680 

375.739 

375.680 

375.680 

x 6 (mm) 

0.828 

0.794 

0.807 

0.846 

0.794 

0.851 

0.848 

0.851 

0.851 

x 7 (mm) 

18.738 

20.000 

20.000 

19.727 

20.000 

20.000 

20.000 

20.000 

20.000 

x 8 (mm) 

4.800 

4.800 

5.000 

4.800 

4.472 

5.000 

5.000 

5.000 

5.000 

x 9 (mm) 

16.000 

16.000 

10.000 

16.000 

10.000 

10.000 

10.000 

10.000 

10.000 


The optimum values of design variables are shown in 
table 5 for each objective function including multiobjective 
functions. 

As shown in table 5, the optimum values of design 
variables are different from one single objective function to 
another. Otherwise, mechanical error f 2 generally induces 
highest values of optimum design variables while their lowest 
values are obtained with manufacturing cost f 5 . It is also 
shown that multiobjective functions f 6 , f 7 , f 8 and f 9 
optimization allow generally to obtain intermediate values of 
these extreme design variables. Therefore, strength reliability 
and mainly manufacturing cost have a considerable effect on 
the optimal design of the middle finger’s mechanism. In fact, 
results obtained from multiobjective functions f 6 , f 8 and f 9 are 
the same and highlight the manufacturing cost role. 

Fig. 10 give optimum values of each objective function 
used in each of optimization problems related to single 
objective and multiobjective functions. 


• From one multiobjective optimization to another, optimum 
values of each single function remain the same whereas 
the fact that optimum values of design variables related to 
optimization function f 6 are different from those of 
optimization functions f 7 , f 8 and f 9 . 

B. Discussion 

In the previous studies [1,2], optimal design variables 
obtained from one multiobjective optimization to another are 
different. In the present work, the global manufacturing cost, 
included in any combination of the simultaneous 
optimization, enabled to obtain the same optimum variables 
design. 

Otherwise, the optimal global manufacturing cost 
obtained in this study is about 495 US Dollars for the middle 
finger mechanism versus 419 US Dollars proposed by 
Ventimiglia [16] which includes complete middle finger 
assembly. Therefore, it will be necessary to evaluate the really 
cost of our complete middle finger assembly. 




■ Single function 

■ f6 (f 1/2/3) 

■ f7 (fl/2/3/4) 

■ f8 (fl/2/3/5) 


Figure 10: Optimum values of objective functions 


It can be shown that: 

• Optimum values of single objective function optimization 
are smaller than those obtained by the means of 
multiobjective function optimization, except those 
related to reliability and manufacturing cost which are a 
little different; 


V. CONCLUSION 

The aim of this study was to introduce into the 
optimization process of the crossed four-bar link mechanism 
of hand prosthesis, two new performance criteria, namely the 
reliability strength of the whole structure and the 
manufacturing cost of the mechanism. These criteria were 
added to those related to the mean quadratic error of the 
bending angle of the second joint of the middle finger, the 
mechanical error due to dimensions tolerances and the 
maximum driving torque applied to the driving shaft. This 
study permits us to confirm that, through the five different 
optimal mechanisms obtained, each criterion used here 
generates a particular design, translating its importance in the 
optimization process. 

Otherwise, because of the influence of each of these 
criteria on the mechanism design, their simultaneously 
optimization showed that the global manufacturing cost, 
included in any combination of the multiobjective 
optimization, enabled to obtain the same optimum variables 
design. Finally, the strength reliability of the driving is always 
equal to the unit. That translated the fact that, the constraints 
induced inside each bar, remain less than the maximum 
material constraints. 

References [17], [18] and [19], related to physic-mechanical, 
chemical and mineralogical characterization of Materials clay 
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of the Far North Region of Cameroon, will be helpful to 
produce ceramic material in order to obtain moulded 
prosthetic fingers in which the manufactured mechanisms will 
be inserted. 
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